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This report describes an observation of mixing-induced CP violation and a measurement of the 
CP violation parameter, sin20i, with the Belle detector at the KEKB asymmetric e + e~ collider. 
Using a data sample of 29.1 fb _1 recorded on the T(45) resonance that contains 31.3 million BB 
pairs, we reconstruct decays of neutral B mesons to the following CP eigenstates: J/ipKg, ip(2S)Kg, 
XciKg, rj c Kg, J/ipK^ and J/ipK . The flavor of the accompanying B meson is identified by 
combining information from primary and secondary leptons, K mesons, A baryons, slow and fast 
pions. The proper-time interval between the two B meson decays is determined from the distance 
between the two decay vertices measured with a silicon vertex detector. 

The result sin20i = 0.99 ± 0.14(stat) ± 0.06(syst) is obtained by applying a maximum likelihood 
fit to the 1137 candidate events. We conclude that there is large CP violation in the neutral B 
meson system. A zero value for sin2^i is ruled out by more than six standard deviations. 



PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 
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I. INTRODUCTION 

The phenomenon of CP violation is one of the major 
unresolved issues in our understanding of particle physics 
today. In 1973, Kobayashi and Maskawa (KM) proposed 
a model where CP violation is accommodated as an ir- 
reducible complex phase in the weak-interaction quark 
mixing matrix ffl , which is defined as 



V u d Vus V u b 

V c d V cs V c b 
V td V u V tb 



(1) 



where the nontrivial complex phases are conventionally 
assigned to the furthest off-diagonal elements V u b and 
Vtd- Unitarity of this CKM matrix (Cabibbo-Kobayashi- 
Maskawa matrix) implies that T,iVijV* k — Sjk, which 
gives the following relation involving V u b and Vtd- 



v ud v: b + v cd v* b + v td v; b = o. 



(2) 



This expression can be visualized as a closed triangle in 
the complex plane as shown in Fig. |} The three interior 
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VcdVcb 

FIG. 1. The unitarity triangle relevant to B decays. Angles 
of the triangle arise from the KM phase. 



angles of the unitarity triangle originate from the non- 
vanishing CP-violating phase (the KM phase) and are 
defined as M: 



b 1 = ir- BXg{-V td V* h f-V cd V; b ) 
b 2 = axg(Vt d V t t/-V ud V: b ), 
b 3 = axg(V ud V: b /-V cd V: b ). 



(3) 



The bold ansatz of Kobayashi and Maskawa required the 
existence of six quarks at a time when only the u, d and 
s quarks were known. The subsequent discoveries of the 
c, b and t quarks as well as the consistency with the 
CP violation observed in the neutral kaon system led to 
the incorporation of the KM mechanism as an essential 
component of the Standard Model (SM) . 

In 1980, Sanda, Bigi and Carter pointed out that 
the KM model contained the possibility of sizable CP- 
violating asymmetries in certain neutral B meson de- 
cays S. The subsequent observation of a long b quark 



lifetime 0] and a large mixing in the neutral B meson 
system [pj indicated that it would be feasible to measure 
CP violation in B meson decays at an asymmetric e + e~ 
collider at the T(AS) energy. 

Until recently the only observation of CP violation was 
in the neutral kaon system, where the interpretation of 
results is complicated due to large corrections from the 
strong interaction. By contrast, these corrections are ab- 
sent or very small for the aforementioned CP violation 
in the neutral B meson system. Thus its measurement 
can be used to over-constrain and test the consistency of 
the SM. 

A pair of neutral B mesons created in the decay 
T(4S) ->■ B°B° is in a state with C = -1 at the time of 
production (t = 0), where C denotes the charge conjuga- 
tion. Although oscillation then starts, the state preserves 
the C odd configuration and is not allowed to be B°B° 
or B Q B°. The time evolution of the pair is given by 



!*(*)) 



- t/T B a Y$>(t = o)), 



(4) 



m = 0)) = ^ [\B (k)B°(-k)) \B Q (-k)B°(k)} 



where k and — k are the B mesons' momenta in the T(4S) 
rest frame. This coherence is preserved until one B meson 
decays. Hence, if we can determine the flavor and the 
decay time itag of one of the B mesons decaying into 
a final state /tag, we are able to determine the time- 
dependent decay amplitude of the other B at any time t 
as a function of the time difference t — t tag . We consider 
the case where the other B meson decays at t — tcp to 
a CP eigenstate, fcp- When the decay is dominated by 
a single transition amplitude, the following formulae for 
the decay rates hold to a good approximation (3): 

R(B° f C p; At) = e-l Ai l/ T BO/ 2TBO 

x [1 + £f sin.2(f>cp siix(Am d At)], 
R(B° -> fcp; At) = e -l At l/^72r B o (5) 

x [1 — £f sia2cj)cp sin(AmdAi)], 

and the time-dependent CP- violating asymmetry is p] 



A(At) = — W 
v ' ulna 



fcp; At) - R(B° f C p; At) 



R(B° -> f C p; At) + R{B° 
= —£f sin 24>cp sm(AnidAt), 



fcp; At) 



(6) 



where £/ is the CP eigenvalue of fcp, Arnd is the mass 
difference between the two B° meson mass eigenstates 
and At = tcp — ttag- Because the asymmetry, A(At), 
vanishes in the time-integrated rate, it is very important 
to measure the time dependence. 

The angle 4>cp is directly related to the interior an- 
gles of the unitarity triangle, and is the phase difference 
between two interfering amplitudes, one for B°(B ) — > 
fcp and the other for the mixing process B°(B°) — > 
B°(B Q ) — > fcp- The quantity <\>cp is equal to <fii if 
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fcp — J/ipKs or an Y other CP eigenstate that arises 
from a b(b) — > ccs(s) transition. The hadronic uncer- 
tainty in this case is negligibly small because the ampli- 
tude of the b — > s flavor-changing transition with asso- 
ciated cc production is not only small but has the same 
weak phase. 

As is described in the next section, the KEKB e + e~ 
collider produces the T(4S) with a Lorentz boost of 
/?7 = 0.425. Since the B° and B° are nearly at rest 
in the T(45*) center of mass system (cms), At can be de- 
termined from the displacement between the two B decay 
vertices — i.e. 



At ~ (z C p - Ztag)//?7 C = Az//?7C, 



(7) 



where the z axis is defined to be anti-parallel to the 
positron beam direction. 

Following initial experimental studies || P] , the BaBar 
fiof and Belle jll| Collaborations recently reported the 
first clear observations of CP violation in the neutral B 
meson system. In this paper, we describe the details of 
the measurement of sin 2</>i with the Belle detector at the 
KEKB asymmetric e + e~ collider with the same 29.1 fb _1 
data sample reported in Ref. |fll]| . In the next section 
we describe the KEKB collider and the Belle detector. 
The measurement of sin 2(j)\ requires the reconstruction 
of B° — > fcp decays (denoted by Bqp), the determi- 
nation of the 6-flavor of the accompanying (tagging) B 
meson, the measurement of At, and a fit of the expected 
At distribution to the measured distribution using a like- 
lihood method. The selection and tagging procedures are 
described in Sections III and IV. After introducing the 



methods to extract sin2^>i from the At distributions in 
Section [v|, we present the results of the fit and discuss 
the interpretation of the CP violation in Section [Vl]. We 
summarize the results in Section VII. 



II. EXPERIMENTAL APPARATUS 

KEKB jl2| is an asymmetric e + e~ collider 3 km in cir- 
cumference, which consists of 8 GeV e~ and 3.5 GeV e + 
storage rings and an injection linear accelerator. It has a 
single interaction point (IP) where the e + and e~ collide 
with a crossing angle of 22 mrad. The data used in this 
analysis were taken between January 2000 and July 2001. 
The collider was operated during this period with a peak 
beam current of 930 mA(e + ) and 780 mA(e~), giving 
a peak luminosity of 4.5 xlO 33 cm _2 s _1 . Due to the en- 
ergy asymmetry, the T(4S') resonance and its daughter B 
mesons are produced with a Lorentz boost of (3j =0.425. 
On average, the B mesons decay approximately 200 /im 
from the T(4S') production point. 

The Belle detector |l3| is a general-purpose large solid 
angle magnetic spectrometer surrounding the interaction 
point. It consists of a barrel, forward and rear compo- 
nents. It is placed in such a way that the axis of the 



detector solenoid is parallel to the z axis. In this way we 
minimize the Lorentz force on the low energy positron 
beam. 

Precision tracking and vertex measurements are pro- 
vided by a central drift chamber (CDC) |l4|] and a sili- 
con vertex detector (SVD) Q. The CDC is a small-cell 
cylindrical drift chamber with 50 layers of anode wires 
including 18 layers of stereo wires. A \cw-Z gas mixture 
(He (50%) and C2H 6 (50%)) is used to minimize multi- 
ple Coulomb scattering to ensure a good momentum res- 
olution, especially for low momentum particles. It pro- 
vides three-dimensional trajectories of charged particles 
in the polar angle region 17° < < 150° in the labora- 
tory frame, where 9 is measured with respect to the z 
axis. The SVD consists of three layers of double-sided 
silicon strip detectors arranged in a barrel and covers 
86% of the solid angle. The three layers at radii of 3.0, 
4.5 and 6.0 cm surround the beam-pipe, a double-wall 
beryllium cylinder of 2.3 cm radius and 1 mm thickness. 
The strip pitches are 84 /jm for the measurement of z co- 
ordinate and 25 fiia for the measurement of azimuthal 
angle <fr. The impact parameter resolution for recon- 
structed tracks is measured as a function of the track 
momentum p (measured in GeV/c) to be a xy — [19 © 
50/(p/3sin 3/2 9)} fim and a z = [36 42/(p/3sin 5/2 9)] /an. 
The momentum resolution of the combined tracking sys- 
tem is a Pt /pt — (0.30//?©0.19p t )%, where p t is the trans- 
verse momentum in GeV/ c. 

The identification of charged pions and kaons uses 
three detector systems: the CDC measurements of 
dE/dx, a set of time-of- flight counters (TOF) [pl| and 
a set of aerog el Cherenkov counters (ACC) [07^ The 
CDC measures energy loss for charged particles with a 
resolution of a(dE/dx) = 6.9% for minimum- ionizing pi- 
ons. The TOF consists of 128 plastic scintillators viewed 
on both ends by fine-mesh photo-multipliers that oper- 
ate stably in the 1.5 T magnetic field. Their time res- 
olution is 95 ps (rms) for minimum-ionizing particles, 
providing three standard deviation (3a) K^/n^ separa- 
tion below 1.0 GeV/c, and 2a up to 1.5 GeV/c. The 
ACC consists of 1188 aerogel blocks with refractive in- 
dices between 1.01 and 1.03 depending on the polar an- 
gle. Fine-mesh photo-multipliers detect the Cherenkov 
light. The effective number of photoelectrons is approx- 
imately 6 for (3 = 1 particles. Using this information, 
P(K/tt) = Prob(K)/(Prob(K) + Prob(n)), the proba- 
bility for a particle to be a meson, is calculated. A 
selection with P(K/ir) > 0.6 retains about 90% of the 
charged kaons with a charged pion misidentification rate 
of about 6%. 

Photons and other neutrals are reconstructed in a 
CsI(Tl) crystal calorimeter (ECL) §J| consisting of 8736 
crystal blocks, 16.1 radiation lengths (Xq) thick. Their 
energy resolution is 1.8% for photons above 3 GeV. 
The ECL covers the same angular region as the CDC. 
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Electron identification in Belle is based on a combina- 
tion of dE/dx measurements in the CDC, the response 
of the ACC, the position and the shape of the elec- 
tromagnetic shower, as well as the ratio of the cluster 
energy to the particle momentum [ pTsj ] . The electron 
identification efficiency is determined from two-photon 
e+e~ — * e+e~e+e~ processes to be more than 90% for 
p > 1.0 GeV/c. The hadron misidentification prob- 
ability, determined using tagged pions from inclusive 
Kg — ► 7r + 7r~ decays, is below 0.5%. 

All the detectors mentioned above are inside a super- 
conducting solenoid of 1.7 m radius that generates a 1.5 T 
magnetic field. The outermost spectrometer subsystem 
is a K® and muon detector (KLM) [^o|, which consists of 
14 layers of iron absorber (4.7 cm thick) alternating with 
resistive plate counters (RPC). The KLM system cov- 
ers polar angles between 20 and 155 degrees. The overall 
muon identification efficiency, determined by using a two- 
photon process e + e~ — > e + e~ /i + fi~ and simulated muons 
embedded in BB candidate events, is greater than 90% 
for tracks with p > 1 GeV/c detected in the CDC. The 
corresponding pion misidentification probability, deter- 
mined using K s — > 7r + 7r~ decays, is less than 2%. 

In our analysis, Monte Carlo (MC) events are gen- 
erated using the QQ event generator |2l| and the re- 
sponse of the Belle detector is precisely simulated by a 
GEANT3-based program |22| . The simulated events are 
then reconstructed and analyzed with the same proce- 
dure as is used for the real data. 

III. RECONSTRUCTION OF B° DECAYS 

We use a 29.1 fb _1 data sample, which contains 31.3 
million BB pairs, accumulated at the T(45 l ) resonance 
between January 2000 and July 2001. The entire data 
sample has been analyzed and reconstructed with the 
same procedure. 

We reconstruct B° decays to the following CP eigen- 
states H): J/ipK%, ip(2S)K%, XaK s and r) c K° s having 
£f = — 1; and J/tpK^ having £f = +1. We also use the 
decay B° -> J/tpK* a , K*° -> K s tt°, which is a mixture 
of even and odd CP eigenstates. The selection of these 
Bqp candidates is described in the following sections. 

A. BB event pre-selection 

To select generic BB candidates, we require at least 
three tracks that satisfy yjx 2 +y 2 < 2.0 cm, \z\ < 
4.0 cm, and pt > 0.1 GeV/c, where x, y, z represent 
the point of closest approach of the track to the beam 
axis, and pt is the momentum of the track projected onto 
the xy-plane. We also require that more than one neutral 
cluster is observed and have energy greater than 0.1 GeV. 



The sum of all cluster energies, boosted back to the cms 
assuming each cluster is generated by a massless particle, 
is required to be between 10% and 80% of the total cms 
energy. The total visible energy in the cms, -E^is" 5 ! * s com- 
puted from the selected tracks, assuming they are pions, 
and the calorimeter clusters that are not associated with 
the tracks. We require that E^ s is greater than 20% of 
the total cms energy. The absolute value of the z com- 
ponent of the cms momentum is required to be less than 
50% of the cms energy. The event vertex reconstructed 
from the selected tracks must be within 1.5 cm and 3.5 cm 
of the interaction region in the directions perpendicular 
and parallel to the z axis, respectively. Monte Carlo sim- 
ulation shows that the selection criteria described above 
retain more than 99% of BB events and J/tp inclusive 
events. 

To suppress continuum background, which consists of 
e + e~ — * qq where q is u, d, s or c quark, we also require 
i?2 = H2/H0 < 0.5, where H2 and Ho are the second and 
zeroth Fox- Wolfram moments p3] . 

B. B° — > Charmonium Kg(K*°) reconstruction 

The candidate J/tp and tp{2S) mesons are recon- 
structed using their decays to lepton pairs, i.e. J/tp — > 
and e + e~. The ip(2S) meson is also reconstructed 
via its J/ipTr + ir~ decay, the \ci meson via its J/tpj decay, 
and the r\ c meson via its K + K~ir° and Kg(Tr + ir~)K ± ir T 
decays. 

For J/tp and tp(2S) — > £ + £~ decays, we use oppo- 
sitely charged track pairs where both tracks are posi- 
tively identified as leptons. For the B° — > J /tpK° s (n + ir~) 
mode, which has the smallest background fraction among 
the CP eigenstates that are used, the requirement for 
one of the tracks is relaxed to improve the efficiency: 
a track with an ECL energy deposit consistent with a 
minimum ionizing particle is accepted as a muon and a 
track that satisfies either the dE/dx or the ECL shower 
energy requirements as an electron. In order to re- 
move either badly measured tracks or tracks that do 
not come from the interaction region, we require \dz\ < 
5 cm for both lepton tracks. In order to account par- 
tially for final-state radiation and bremsstrahlung, the 
invariant mass calculation of the e + e~ pairs is corrected 
by adding photons found within 50 mrad of the e + 
or e~ direction. Nevertheless, a radiative tail remains 
and we use an asymmetric invariant mass requirement 
-150 < M e+e - - Mj/^(2S)) < 36 MeV/c 2 . Since 
the radiative tail is smaller, we select —60 < 

M„ +JU - - M J/MM2S)) < 36MeV/c 2 ||. Events with 
a candidate J/ip — > £ + £~ decay are accepted if the J/tp 
momentum in the cms is below 2 GeV/c. Figure ^ shows 
the invariant mass distributions for J/tp — > fi + and 
J/tp — > e + e _ with the selection criteria applied for the 
J/tpKg mode. 
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FIG. 2. The invariant mass distributions for (a) 
J/tp — > an d (b) J/tp — > e + e~ in the reconstruction 

of the J/tp Kg mode where the selection for one of the tracks 
is relaxed (details are explained in the text). 



To reconstruct ip(2S) — > J/i/j7r + n~ decays, we se- 
lect n + ir~ pairs with an invariant mass greater than 
400 MeV/c 2 . This requirement is based on the mea- 
sured mass distribution J2(| and improves the signal-to- 
background ratio. The tp(2S) candidates are then se- 
lected requiring the mass difference, Mg+t-^-K- —Mg+i- , 
to be between 0.58 GeV/c 2 and 0.60 GeV/c 2 . This cor- 
responds to a ±3<7 requirement where a is the resolution 
on the mass difference . 

The Xci — * J/^l candidates are selected by requir- 
ing the mass difference, M^+^- 7 — M e + e ~ , to be between 
0.385 GeV/c 2 and 0.4305 GeV/c 2 . We veto photon can- 
didates that form a good ir° candidate with any other 
photon candidate of energy greater than 60 MeV in the 
event. A good ir° candidate is defined by an invariant 
mass within —28 to +17 MeV/ c 2 of the nominal ir° mass, 
and by a x 2 of less than 10 after a mass-constrained kine- 
matic fit 

For Kg — * ir + n~ reconstruction, we select oppositely 
charged track pairs that satisfy the following require- 
ments: (1) when both pions have associated SVD hits, 
the distance of closest approach of both pion tracks in 
the z direction should be smaller than 1 cm; (2) when 
only one of the two pions has associated SVD hits, the 
distance of closest approach of both the pion tracks to 
the nominal interaction point in the x-y plane should 
be larger than 0.25 mm; (3) when neither pion has an 
associated SVD hit, the <p coordinate of the 7r + 7r~ ver- 



tex and the <p direction of the 7r + 7r~ candidate's three 
momentum vector should agree within 0.1 radian. The 
invariant mass of the candidate 7r + 7r~ pair is required to 
be between 482 and 514 MeV/c 2 , which retains 99.7% of 
the Kg candidates. 

For the XciKg and rj c Kg modes, more stringent track 
selection criteria are applied in Kg — > ir + ir~ reconstruc- 
tion to reduce the background: (1) the flight length in 
the r-(p plane should be greater than 1 mm (2 mm for 
77c — > K + K~tt°); (2) a mismatch in the z direction at 
the Kg vertex point for two 7r ± tracks should be less 
than 2.5 cm (10 cm for r] c Kg); (3) the angle in the r-(p 
plane between the Kg momentum vector and the direc- 
tion defined by the K s and J/tp (or rj c ) decay vertices 
should be less than 0.2 (0.1 for r] c Kg) radian; and (4) for 
the XciKg selection we also require that the distance of 
closest approach of the J/ip vertex in the radial direction 
for each tt track should be greater than 0.25 mm. 

To reconstruct Kg — > ir°ir° candidates, we first se- 
lect photons that have an energy of at least 20 MeV. 
For 7T° — * 77 candidates, we require that the invari- 
ant mass of the two photons be between 80 MeV/c 2 
and 150 MeV/c 2 and the momentum of 7r° be greater 
than 100 MeV/c. Initially, the 7r° decay vertex is as- 
sumed to be at the nominal interaction point. To select 
Kg — > 7r°7r° candidates, we find the best decay vertex 
where the invariant masses of two ir° candidates are the 
most consistent with the nominal ir° mass. To this end, 
first the Kg flight direction is measured from the sum of 
the momenta of the four photons, then we calculate the 
X 2 of the mass constrained fit for each 7r°, varying the 
decay vertex along the Kg direction through the IP. We 
choose the vertex point that minimizes the sum of the x 2 
for the two ir° candidates. We require that the distance 
between the IP and the reconstructed Kg decay vertex 
be larger than —20 cm where the positive direction is 
defined by the Kg momentum. To reduce the combina- 
torial background, the x 2 f° r each tt° meson at the Kg 
decay vertex point determined by this method is also re- 
quired to be less than 10. Using the calculated Kg decay 
vertex, we finally require that the invariant mass of Kg 
candidate lie between 470 MeV/c 2 and 520 MeV/c 2 . 

For J /ipK*°(KgTT°) decays, we use K s Tr° combina- 
tions that have an invariant mass within 75 MeV/c 2 of 
the nominal K*° mass. Here, the 7r° candidate is re- 
constructed from photons with an energy greater than 
40 MeV and the two photon invariant mass is required 
to be between 125 and 145 MeV/c 2 . We reduce the 
background from low-momentum ir° mesons by requiring 
cos8k* < 0.8, where Ok- is the angle between the K*° 
flight direction and the Kg momentum vector calculated 
in the K*° rest frame. 

We reconstruct B° — ► rj c Kg candidates in two r\ c de- 
cay modes: T) c — > KgK ± ir T and r/ c — > K + K~ir°. We 
require charged kaons to be positively identified using 
CDC dE/dx measurements and information from the 
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TOF and ACC systems. For the r] c -> K^K^ir? chan- 
nel, we require an invariant mass ranging from 2.935 to 
3.035 GeV/c 2 . In order to suppress the continuum back- 
ground, we require i?2 < 0.45 and | cos # t hr| < 0.85, where 
#thr is the angle between the thrust axis of the B° candi- 
date and that of all remaining charged and neutral par- 
ticles in the event. In the r\ c — ► K + K~ir° mode, we 
reconstruct the it meson from photons having an energy 
larger than 50 (200) MeV in the ECL barrel (end-cap) 
region. The invariant mass of the r\ c — ► K + K~ir° candi- 
date is required to be between 2.890 and 3.040 GeV/c 2 . 
The continuum background is suppressed by requiring 
R 2 < 0.40 and | cos0 t hr| < 0.80. 

For B reconstruction, we calculate the energy dif- 
ference, AE, and the beam-energy constrained mass, 
Afbc- The energy difference is defined as AE = E^ ms — 
and the beam-energy constrained mass Mb c = 
lere E™* m is the cms beam en- 
are the cms energy and momen- 
To improve the momentum 



pcms 
beam 



ergy, and E™ s and p™ 
turn of the B° candidate, 
resolution, a vertex fit and then a mass-constrained fit 
are performed wherever needed. The resulting fitted mo- 
menta are used in the AE and Mb c calculations. A scat- 
ter plot of Mbc and AE for J/4>Kg(Tr + Tr~) candidates 
is shown in Fig. along with the projections onto each 
axis. The B° candidates are selected by requiring 5.270 

< M bc < 5.290 GeV/c 2 (|M bc - M B o\ < 3.5cr) and by 
applying the mode-dependent requirements on AE listed 
in Table H. Figure ^ shows the Mb c distribution after the 
AE selection. 

The background in the signal region is estimated by 
simultaneously fitting the Md C and AE distributions 
with signal and background functions in the region 5.2 

< M bc < 5.3 GeV/c 2 and -0.1 < AE < 0.2 GeV. We 
use a two-dimensional Gaussian for the signal. For the 
background, we use the ARGUS background function J2^] 
for Mbc and a linear function for AE. The details are 
described in Section 



VD 



The number of candidates ob- 
served as well as the estimated background are given in 
Table |. 



C. B° — ► J/xpK" reconstruction 

The reconstruction of B° — > J/xpKl is an experimen- 
tal challenge but is very important because its yield is 
expected to be large. In addition, since this mode is a 
CP-even eigenstate, we should observe a time-dependent 
asymmetry reversed in sign compared to J/tpKg, which 
provides an important experimental consistency check. 

While the detached vertex and invariant mass of the 
Kg provide significant background reduction for J/ipKg, 
the background is larger for J/ipK^ as only the K\ di- 
rection is measured. Since the energy of the K® L is not 
measured, Mb c and AE cannot be used as the final kinc- 
matical variables to identify B° candidates as in other fi- 
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FIG. 3. The scatter plot of AE versus M bc for 
J/ipKg(ir + iv~) candidates. The box represents the signal re- 
gion. The upper left figure is the AE projection with 5.270 
< M bc < 5.290 GeV/c 2 . The lower right figure is the M hc 
projection with \AE\ <0.04 GeV. The enhancement in the 
negative AE region is due to decay modes with additional 
pions, e.g. B -> J/ipK*, K* -> K° s n. 
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FIG. 4. The Mbc distribution for all exclusive decay 
modes combined except J/^K%. The AE cuts are imposed. 
The shaded area is the estimated background. 



nal states. Using the four-momentum of a reconstructed 
J/tp candidate and the K\ flight direction, we calculate 
the momentum of the K® candidate requiring AE = 0. 
We then calculate p™ s which is used for the final selec- 
tion. 

The selection criteria are necessarily tighter than those 
used for the J/ipKg candidates. However, precise deter- 
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TABLE I. Summary of AE signal region, the numbers of 
signal candidates, and expected background. 



Decay 


AE cut (MeV) 


Signal Expected 


mode 

nil i/.ltA) 


Lower ~ Upper 


candidates background 



K% -> 7T+7T- 


-40 - 


- 40 


457 


11.9 


Kg — > 7T n 


-150 - 


- 100 


76 


9.4 


B° -> 4>(2S)K S 










V>(25) -> i+r 


-40 - 


- 40 


39 


1.2 




-40 - 


- 40 


46 


2.1 


B° - Xe i*S 


-40 - 


- 40 


24 


2.4 














-40 - 


- 40 


23 


11.3 


r; c -> K + K~1T° 


-60 - 


- 40 


41 


13.6 


S° -> J/V>iT° 


-50 - 


- 30 


41 


6.7 



mination of the Kl flight direction with the KLM and 
ECL allows us to reconstruct J/ipK^ candidates with 
sufficient efficiency and purity. 

We use tracks which are positively identified as 
electrons (muons) in the identification of J/ip — » 
e + e~ (/i + /i~) candidates. We require the invariant mass 
of the lepton pair to lie in the range 3.05 < M e + e - < 
3.13 GeV/c 2 . The radiative photon correction for elec- 
tron pairs is made in the same way as for the other modes. 
Events are rejected if one of the following decay modes 
are exclusively reconstructed, and satisfy \AE\ < 0.05 
GeV and 5.27 < M hc < 5.29 GeV/c 2 : J/ipK+, J/ipK% 
J/tpK* + (K*+ -> K+n ,K^TT+), and J/ipK*° {K*° -> 
K+7r-,K a s 7T°). 

We select K\ candidates based on the KLM and ECL 
information. There are two classes of K\ candidates that 
we refer to as KLM and ECL candidates. To select the 
KLM candidates, a cluster of KLM hits is formed by 
combining the hits within a 5° opening angle. We require 
hits in two or more KLM layers and calculate the center of 
the KLM cluster. If there is an ECL cluster with energy 
greater than 0.16 GeV within a 15° cone, we relax our 
criteria to allow a cluster with a hit in just one KLM 
layer. In this case the direction of the ECL cluster is 
taken as the K a L direction. If the cluster lies within a 15° 
cone of the extrapolation of a charged track to the first 
layer of the KLM, it is discarded. 

ECL candidates are selected from ECL clusters using 
the following information: the distance between the ECL 
cluster and the closest charged track position; the ECL 
cluster energy; the ratio of energies summed in 3 x 3 
and 5x5 arrays of Csl crystals surrounding the crystal 
at the center of the shower; the ECL shower width and 
the invariant mass of the shower. After a very loose pre- 
selection based on the above five discriminants, we calcu- 
late signal and background likelihood values for each dis- 
criminant based on a J/ip inclusive MC. Taking the prod- 
ucts of the above five likelihoods for each signal and back- 



ground, we form the likelihood ratio C K o / (C K o + £f a kc)- 
This likelihood ratio is required to be greater than 0.5. 

We examine the characteristics of iC£ candidates in 
both the data and the MC. We obtain consistent distri- 
butions for the number of candidates per event, the 
K\ flight directions in the laboratory system, the total 
number of hit KLM layers and the number of hit first- 
layers in the K\ candidates. We investigate the momen- 
tum dependence of the KLM response by using charged 
pions and kaons. The MC simulation reproduces the re- 
sults obtained with data well. We also use e + e~ — > j<p 
followed by — > K ( [Kg, where the exact K\ direction 
and momentum can be obtained by reconstructing 7 and 
Kg — > 7r + 7r~. These studies indicate that the if£ identi- 
fication is well reproduced by the MC with an exception 
of an overall detection efficiency. The K\ detection ef- 
ficiency in the data is found to be lower than the MC 
expectation. This, however, does not cause a difficulty 
in our analysis since we do not rely on the MC de- 
tection efficiency. 

For both KLM and ECL candidates, we require that 
the K\ direction should be within 45° of its expected 
direction calculated from the J/ip candidate momentum 
assuming that the B° candidate was at rest in the cms. 
We also require that no photon from a ir° decay is found 
near the candidate. For this requirement, we select 
7T° candidates satisfying 0.12 < M 77 < 0.15 GeV/c 2 and 
with momentum above 0.8 (1.2) GeV/c for KLM (ECL) 
candidates. 

To reconstruct B° -> J/ipK° L , we first use the KLM 
candidates. If none of the KLM candidates satisfy the se- 
lection criteria below, we use the ECL candidates. Thus, 
the two classes of B Q — > J/ipK^ candidates are mutu- 
ally exclusive. In order to suppress the background, we 
calculate a probability density function (pdf ) for each of 
the following variables and then form a product of the 
pdf's to obtain the combined likelihood: the cms mo- 
mentum of the J/ip, P C JJ^', the angle between the K\ 
candidate and the closest charged track having a mo- 
mentum larger than 0.7 GeV/c; the Fox-Wolfram mo- 
ment ratio R2; cos 9 b where b is the polar angle of the 
reconstructed B in the cms; and the number of charged 
tracks with p t > 0.1 GeV/c, \dr\ < 2 cm, and \dz\ < 
4 cm. In addition to these five variables, two other vari- 
ables are included conditionally to reduce the background 
from B+ -> J/ipK*+, K*+ -> K° l tt + decays. One such 
variable is the cms momentum of the J/iPK ( /tt + system, 
p c ^ as (J/ipK l / j Tr), and the other is the momentum of the 
additional pion. We use p c B ms ( J/ipK^jr) if the invariant 
mass of the and a charged track, with the nominal 
pion mass being assumed, is above 0.85 GeV/c 2 and be- 
low 0.93 GeV/c 2 , and p^ ms ( J/ipKlw) < 0.8 GeV/c |j§. 
The pion momentum is used if the addition of the extra 
pion results in 0.2 < p c ™ (J / ip K° L ir) < 0.45 GeV/c given 
the requirement above on the K^ir^ invariant mass. 

The combined likelihood is calculated for J / ipK\ (sig- 
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nal) and inclusive B — > J/tp decays (background) us- 
ing MC samples. Background events from misidentified 
(fake) J/tp mesons are not included in the likelihood con- 
struction since their contribution is small. We then form 
a likelihood ratio Lj^ K o /{Cj^ K o + £bkg) that is used 
as a discriminant variable. Figure ^| shows the likelihood 
ratio distributions for the data and Monte Carlo candi- 
dates. We require that the likelihood ratio be larger than 



parameter of the fit is an offset in p c g ns , allowing the 
signal shape to shift with respect to the background dis- 
tribution. 

The result of the fit to the p^ 11 * 3 distribution is shown in 
Fig. ^|. By integrating each component obtained by the 
fit in the signal region, we find a total of (346.3 ± 28.8) 
J/tpK 1 / signal events, and a signal purity of 61%. The 
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FIG. 5. The likelihood ratio distributions for (a) J/tpK% 
Monte Carlo, (b) J/tp + X background Monte Carlo, and (c) 
the data in the signal region with Monte Carlo events overlaid. 



0.4 to identify J/tpK ( [ candidates. 

Figure || shows the resultant p™ s distribution for 
events that satisfy all the selection criteria. We define 
the signal region as 0.2 < pf 18 < 0.45 (0.40) GeV/c 
for K\ candidates identified with the KLM (ECL) crite- 
ria p0[, and obtain 569 candidates, out of which 397 are 
KLM candidates and 172 are ECL candidates. 

We extract the J / ipK® signal yield by fitting the p c g ls 
distribution of the data to a sum of four components: (1) 
signal; (2) background with K\; (3) background with- 
out K\\ and (4) combinatorial J/tp mesons. The shapes 
of the first three components are determined from the 
J/tp inclusive MC and look-up tables are used in the fit. 
The normalizations of these three components are treated 
as free parameters in the fit to minimize the effect of 
the aforementioned uncertainty in the K® detection effi- 
ciency in the MC simulation. The combinatorial compo- 
nent is evaluated using events with e/i-pairs that satisfy 
the requirements for J/tp reconstruction. The shape is 
modeled by a second-order polynomial. An additional 
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FIG. 6. The p c ^ a distribution of the J/tpK° L candidates. 
The solid line is a sum of the signal and background. The 
shaded histogram shows the background component only. 



fitting procedure is applied to the KLM and ECL candi- 
dates separately. The \ 2 of the fit to KLM (ECL) can- 
didates is 54.8 (31.4) for ndf = 35. Table || summarizes 
the fit results. 



TABLE II. Results of the fit to the p™ 8 distribution of the 
J/tpK% candidates. Yields are for the signal p™ 5 region. 





KLM 


ECL 


Sum 


Signal 


244.8 ± 25.1 


101.5 ± 14.2 


346.3 ± 28.8 


Background 








with Kl 


118.1 


26.5 


144.6 


without Kl 


27.5 


38.6 


66.1 


Combinatorial 


6.6 


5.4 


12.0 


Total Yield 


397 


172 


569 



D. Control samples of flavor-specific B decays 

The reconstruction of flavor-specific B decays is also a 
key ingredient in this analysis, since such decays are used 
for various purposes: evaluation of the performance of fla- 
vor tagging (Section IV B); extraction of the proper-time 
resolution function (Section |VB|); B lifetime measure- 
ments as a cross check (Section VI B 2j ) ; and the demon- 
stration of null asymmetries in non-CP final states (Sec- 
tion VI B 3 ) . Since a large number of events with high pu- 
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rity are required for these purposes, we use semileptonic 
decays and hadronic decays from b — > cud transitions. 

For the semileptonic decays, we use the decay chains 
B° -> D*-(+v and D*~ -> TP n~ , where D° -> if+TT", 
U° -> iC+Tr-vr and Z>° -> K+tt+tt-tt-. All tracks used 
must have associated SVD hits and radial impact pa- 
rameters \dr\ < 0.2 cm, except for the slow pion from the 
D*~ decay. Candidate D Q decays are selected by requir- 
ing the invariant mass to be consistent with the nominal 
D mass. The invariant mass requirement depends on 
the D° decay mode, varying from +9 to +23 MeV/c 2 
above the D° mass and from —9 to —37 MeV/c 2 be- 
low the D° mass, respectively. For the D*~ reconstruc- 
tion, we also impose a requirement on the mass differ- 
ence, Mdiff, between a D*~ candidate and the corre- 
sponding D° candidate. We require Mdiff to be within 
0.8 to 1.75 MeV/c 2 of the nominal mass difference, de- 
pending on the D° decay mode. In addition, we require 
the D*~ cms momentum to be less than 2.6 GeV/c to 
suppress D*~ mesons from the continuum. The D*~ 
candidates are combined with p + or e + candidates hav- 
ing the opposite charge to the D*~ candidate. Lepton 
candidates must satisfy 1.4 < pf ns < 2.4 GeV/c, where 
p^ ms is the cms momentum of the lepton. We exploit the 
massless character of the neutrino, to calculate M^ iss , 
the effective missing mass squared in the cms defined 
by M 2 iss ee (Ef™> ~ £g55) 2 - \p™f - \v™i\\ and 
a product of the momenta of the B and the D*l sys- 
tem, C = 2|pJ ms | \p^f\. The cosine of the angle be- 
tween Pg ms and p™l is given by —M^ iss /C. In the 
Af^ iss versus C plane, therefore, we select D*tv candi- 
dates inside the triangle shown in Fig. |j| The trian- 



gle is defined by C > 



1-45 yfl (-< > 
1.595 miss' w — 

1.45. 
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The vertex position of the D*lv candidate is obtained 
by first fitting the vertex of the D candidate from its 
daughter tracks, and then fitting the lepton and D can- 
didate to obtain the B vertex. The slow pion track from 
the D* meson is also included in the fit. This is im- 
portant since the re-fitted helix parameters of the slow 
pion are used to recalculate Mdiff and improve its reso- 
lution. The overall signal fraction in these control sam- 
ples is estimated to be 79.4%. The backgrounds consist 
of fake D* mesons (10.4%), incorrect combinations of 
D* with leptons that do not show an angular correla- 
tion (4.1%, called "uncorrelated background"), contin- 
uum events (2.1%) and B^ decays (4%). The fraction 
of fake leptons in the uncorrelated background is esti- 
mated with a Monte Carlo simulation to be 0.5% and is 
neglected. We estimate the fake D* background by us- 
ing events in the D° mass sideband as well as fake D* 
events that are reconstructed with wrong-charge slow pi- 
ons. We evaluate the uncorrelated background fraction 
by counting signal events in a sample wherein we flip the 
candidate lepton momentum vector artificially. In this 
case the number of uncorrelated background in the sig- 
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FIG. 7. The distribution of candidate B° -> D*-£ + v de- 
cays in the plane of C versus M^ iss . The signal region (tri- 
angle) is also shown. 



nal region remains the same level while the signal events 
are rejected We apply the same event selection to 

off-resonance data ( 2.3 fb _1 ) and count the number 
of events in the signal region after subtracting the fake 
D* background. We estimate the continuum background 
fraction by scaling the result with the integrated luminos- 
ity. We fit the C distribution with a range — 10 < C < 1.1 
to estimate the signal and B — > D**£u background frac- 
tions. We use MC information to model the signal and 
B — > D**£v distribution. The C distributions and frac- 
tions of all the other backgrounds are obtained from the 
aforementioned special background samples and are fixed 
in the fit. We treat the background fraction uncertain- 
ties as a source of systematic errors in the determination 
of wrong tag fractions, which will be discussed in Sec- 
tion |IV | . 

For hadronic decays, we use the decay modes B — > 
D~n + , D*~n + and D*~p + , where D~ — > K + tt~it~, 
p + — > 7T + 7r°. We reconstruct D*~ candidates in the same 
modes that were used for the D*~£ + T> mode. For D and 
D*~ candidates we apply mode-dependent requirements 
on the reconstructed D mass (ranging from ±30 to ±60 
MeV/c 2 ) and M dm (ranging from ±3 to ±12 MeV/c 2 ), 
in a similar way as for D*~ l + v mode. We select p + 
candidates by requiring the 7r + 7r° invariant mass to be 
within 150 MeV/c 2 of the nominal p + mass. In order 
to suppress continuum background, we impose mode de- 
pendent cuts on i?2 (upper cut values ranging from 0.5 
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to 1.0) and cos# t hr (upper cut values ranging from 0.92 
to 1.0). The cut values are chosen to maximize the fig- 
ure of merit S/\/S + B for each mode, where S and B 
are numbers of signal and background, respectively. We 
select B° candidates by requiring 5.27 < Mb c < 5.29 
GeV/c 2 and \AE\ < 50 MeV. Background contributions 
are estimated by fitting the Mb c and AE distributions 
in the same way as we do for CP modes. Figure || shows 
the Mb c distribution for the three control sample modes 
combined. The overall signal purity is ~82%. We study 
background components with a MC sample that includes 
both BB and continuum events. We find no significant 
peaking background in the signal region. Therefore we 
use the ARGUS function to model the background in the 
Mbc distribution. A possible deviation due to combinato- 
rial backgrounds may introduce ambiguities in the signal 
fraction determination and the background At model pa- 
rameters of the control sample, resulting in small uncer- 
tainties in the final determination of the CP asymmetry. 
This is considered as a source of systematic uncertainties. 



TABLE III. Summary of the numbers of flavor-specific de- 
cay candidates and their estimated purities. 
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FIG. 8. Mbc distribution for all hadronic decay modes 
(5° -» D--n+, D*--n + , and D*~ p+) with \AE\ < 50 MeV. 
Fit curves for background (dashed) and signal plus back- 
ground events (solid) are superimposed. 



The numbers of candidates and their purities for flavor- 



specific samples are summarized in Table [II 



IV. FLAVOR TAGGING 

After the exclusive reconstruction of a neutral B meson 
decaying into a final state, /cp, all the remaining par- 
ticles should belong to the final state, / tag , of the decay 
of the other B meson. To observe time-dependent CP 



Decay mode 


Candidates 


Purity 


B u -> D*-e+v 


16101 


79.4 % 


B u D~n+ 


2241 


85.5 % 


B" -> D*~n+ 


2126 


87.2 % 


B a D*-p+ 


1620 


72.2 % 



violation, we need to ascertain whether / tag is from a B° 
or B°. This determination is called "flavor tagging." The 
simplest and most reliable method for flavor tagging uses 
the charge of high- momentum leptons in semileptonic de- 
cays, i.e. B° -> X£+v and B° -> Xt~V. The charges 
of final-state kaons can also be used since the decays 
B° -> K+X (with b -> c -> s) and B° -> K~X (with 
b —* c — > s) dominate. In addition to these two lead- 
ing discriminants, our algorithm includes other categories 
of tracks whose charges depend on the b quark's flavor: 
lower momentum leptons from c — > sl + i>; A baryons from 
the cascade decay b — ► c — * s; high-momentum pions 



that originate from decays like B° — > D^*'(n + , p + , af, 
etc.); and slow pions from D*~ — ► D a n~. All these 
inputs are combined, taking their correlations into ac- 
count, in a way that maximizes the flavor tagging per- 
formance. The performance is characterized by two pa- 
rameters, e and w. The parameter e is the raw tagging 
efficiency, while w is the probability that the flavor tag- 
ging is wrong (wrong tag fraction). A non-zero value of 
w results in a dilution of the true asymmetry. For ex- 
ample, if the true numbers of reconstructed B° and B° 
are n B o and ro-g > the corresponding asymmetry is A. = 
(n bo — n-g )/{ngo + n-g ) . With realistic flavor tagging, 
the observed numbers are Nga = e((l — w)riBa + wn-g ) 
for B° , iV-g = e((l — w)n-g + umgo) for B°, and the 
observed asymmetry becomes (1 — 2w)A. Since the sta- 
tistical error of the measured asymmetry is proportional 
to e -1 / 2 , the number of events required to observe the 
asymmetry for a certain statistical significance is propor- 
tional to e ff = e(l — 2u;) 2 , which is called the "effective 
efficiency." Note that an imperfect knowledge of w shifts 
the central value of the measurement and thus represents 
a potential source of systematic error. 

In light of the above, our tagging algorithm has been 
designed to maximize e c ff- Moreover, since w directly 
affects the central value of our result, we have developed 
an approach wherein it can be determined from the data. 

In our approach, we use two parameters, q and r, to 
represent the tagging information. The first, q, corre- 
sponds to the sign of the b quark charge where q = +1 
for b and hence B°, and q = — 1 for b and B°. The 
parameter r is an event-by-event flavor-tagging dilution 
factor that ranges from r = for no flavor discrimination 
to r = 1 for unambiguous flavor assignment. The values 
of q and r are determined for each event from a look-up 
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table. Each entry of the table is prepared using a large 
statistics MC sample and is given by 



N{B a ) - N(B°) 
N(B°) + N(B°) ' 



(8) 



where N(B°) and N(B°) are the numbers of B° and B° 
in the MC sample. 

In this analysis, we sort flavor-tagged events into six 
bins in r. For each r bin, we empirically determine w 
directly fr om da ta by using control samples, as described 
in Section IV B below. 



A. Flavor Tagging Method 

Flavor tagging proceeds in two stages. In the first 
stage, the flavor tagging information (q and r) provided 
by each track in the event is calculated. In the second, 
the track-level results are combined to determine event- 
level values for q and r. 

Tracks are sorted into four categories, namely those 
that resemble leptons, kaons, A baryons, and slow pions. 
For each category, we consider several tagging discrimi- 
nants, such as track momentum and particle identifica- 
tion information. The value of q and r for each track 
is assigned based on MC-generated look-up tables that 
take the tagging discriminants as input. 

In the second stage, the results from the four track 
categories are combined to determine the values of q and 
r for each event. Again a look-up table is prepared to 
provide q ■ r. 

Figure || shows a schematic diagram of the flavor tag- 
ging method. The event-level parameter r should satisfy 
r ~ 1 — 2w where we measure w from control samples. 
Using this MC-determined dilution factor r as a measure 
of the tagging quality is a straightforward and powerful 
way of taking into account correlations among various 
discriminants. Using two stages, we keep the look-up ta- 
bles small enough to provide sufficient MC statistics for 
each entry. In the following, we provide additional details 
about each stage of the flavor tagging. 



1. Track-level flavor tagging 

We select tracks that do not belong to Bqp and that 
satisfy \dr\ < 2 cm and \dz\ < 10 cm. Tracks that are 
part of a Kg candidate are not used. Each selected 
tag-side track is examined and assigned to one of the 
four track categories. Tracks in the lepton category are 
subdivided into categories for electron-like and muon- 
like tracks. If the cms momentum, pf ns , of a track is 
larger than 0.4 GeV/c and the ratio of its electron and 
kaon likelihoods is larger than 0.8, the track is assigned 
to the electron-like category. If a track has p c ™ s larger 
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FIG. 9. A schematic diagram of the two-stage flavor tag- 
ging. See the text on the definition of the parameters "5" and 



than 0.8 GeV/c and the ratio of its muon and kaon like- 
lihoods is larger than 0.95, it is passed to the muon- 
like category. The likelihood is calculated by combining 
the ACC, TOF, dE/dx, and ECL or KLM information. 
In the lepton category, leptons from semileptonic B de- 
cays yield the largest effective efficiency. Leptons from 
B — > D cascade decays and high-momentum pions from 
B° — > D^*^n + X also make a small contribution to this 
category. We choose the following six discriminants: the 
track charge; the magnitude of the momentum in the 
cms, p™ s ; the polar angle in the laboratory frame, #i a t,; 
the recoil mass, M reC oii, calculated using all the tag side 
tracks except the lepton candidate; the magnitude of the 
missing momentum in the cms, -P^si anc ^ the lepton-ID 
quality value. 

The track charge directly provides the 6-flavor q. The 
lepton-ID quality distinguishes leptons from pions. Its 
performance is reinforced by variables p^ ms and Ouh, 
which have distributions that are different for leptons and 



pions. The variables p^ ms , M reco ;i and discriminate 
between high momentum leptons from semileptonic B de- 
cays and intermediate momentum leptons from B — > D 
cascade decays where the D decays semi-leptonically. 

If a track cannot be positively identified as a kaon and 
its momentum is less than 0.25 GeV/c, it is assigned to 
the slow-pion category, since low-momentum pions often 
come from charged D* — > Dir decays. Here the discrimi- 
nant variables are: the track charge; the momentum and 
polar angle in the laboratory frame, pi a b and #i a b; the 
ratio of the electron to tt probability from dE/dx, and 
cosathr, the cosine of the angle between the slow pion 
candidate and the thrust axis of the tag-side particles in 
the cms. The main background in this category comes 
from other (i.e. non-Z?* daughter) low momentum pions, 
electrons from photon conversions and ir° Dalitz decays. 
To separate slow pions from those electrons we use dE/dx 
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for this class. Since the direction of the slow pion from a 
D* decay is approximately parallel to the D* direction, 
it is also almost parallel to the thrust axis. The vari- 
ables cosathr, Piab and #i a b, thus, help to identify the 
slow pions originating from D* decays. 

If a track forms a A candidate with another track, it 
is assigned to the A category. In that category the dis- 
criminant variables are: the flavor (A or A) ; the invariant 
mass of the reconstructed A candidate; the angle differ- 
ence between the A momentum vector and the direction 
of the A vertex point from the nominal IP; the mismatch 
in the z direction of the two tracks at the A vertex point; 
and the proton-ID quality value. 

If a track does not fall in any of the categories described 
above, and is not positively identified as a proton, it is 
classified as a kaon. The kaon category is subdivided into 
two parts, one for events with Kg decays, and the other 
for events without Kg's. Separate treatment is necessary 
as events with K$ have a larger wrong tag fraction be- 
cause of their additional strange quark content. We use 
the track charge, p cms , #i a b and the probability ratio of 
kaon to pion as the tagging discriminants. The charge 
of kaons is the most important discriminant. The other 
three variables help separate kaons from pions. 

Although the discriminating power of high-momentum 
pions is weaker than that of charged kaons, they do pro- 
vide some tagging information and are therefore included 
in the kaon category. Approximately half of the pions 
with p cms > 1.0 GcV/c are included in the kaon category, 
while the other half falls into the lepton class, mostly in 
the muon-like category. 

As is described in Section IV E| , we have developed 
a method wherein the wrong tag fractions in our flavor 
tagging method are evaluated from data. Thus possible 
discrepancies between data and MC in the distributions 
of discriminant variables do not affect our sin 20i mea- 
surement, although they might result in the degradation 
in the effective efficiency Nevertheless, we have made 
detailed comparisons between data and MC, which are 
described elsewhere |||, and obtain consistent distribu- 
tions. 



2. Event-level flavor tagging 

The event-level flavor tagging combines the results 
from each of the track categories to determine an overall 
q and r. For the lepton and slow-pion track categories, 
we take the ^-flavor assignment from the track with the 
highest r-value in each category. For the kaon and A cat- 
egories, a combined ^-flavor output is calculated as the 
product of likelihood values for all tracks: 



where the subscript i runs over all tracks in the kaon and 
A categories. The product likelihood is designed to use 
the information from the sum of the strangeness, which 
provides better flavor-tagging performance than simply 
choosing the best candidate. 

Using the three aforementioned track-level q ■ r val- 
ues, the event-level q and r values are obtained from a 
look-up table that is prepared with a MC sample that is 
independent of the sample used to obtain q ■ r values in 
the track categories. The probability that we can assign 
a non-zero value for r is 99.6% in MC; i.e. almost all the 
reconstructed candidates can be used to extract sin 2(j>\ . 

We specify the following six regions : < r < 0.25, 
0.25 < r < 0.5, 0.5 < r < 0.625, 0.625 < r < 0.75, 
0.75 < r < 0.875 and 0.875 < r < 1. For each region we 
obtain the wrong tag fraction, wi, where I is the region 
ID (I = 1,..., 6), using hadronic and semileptonic con- 
trol samples, which is described in the next section. In 
this way, the analysis is insulated from systematic differ- 
ences between the MC simulation and the data due to 
imperfections in the modeling of the detector response, 
decay branching fractions, and fragmentation in our MC 
simulation. 

As a validation, we compare the distribution of q ■ r in 
the D*iv control sample with the MC expectation. As 
shown in Fig. [l0|, the data and MC are in good agree- 
ment. 



1200 
1000 
800 
600 
400 
200 





1200 
1000 
800 
600 
400 
200 




(a) 



B u 



1 -0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 

q-r 



(b) 



B u 



i ... i ... i ... i ... i ... i . . 



1 -0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 

q-r 



(q-r) 



if/A 



nji + (g-^]-nji-(g-od 

IL[i + («•»•)<] +IL[i- («• »•)<]• 



(9) 



FIG. 10. The q ■ r distribution for (a) B° and (b) B° can- 
didates. In each figure the closed (open) points with errors 
show the B° — > D*~ l + v (hadronic B° decays) data with the 
background subtracted, while the histogram is the MC pre- 
diction. 
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B. Flavor Tagging Performance 



The flavor tagging performance is evaluated by replac- 
ing the CP-cigenstate side of the event with a flavor- 
specific decay and tagging the 6-flavor for the other side 
using the method described above. We use the semilep- 
tonic decay B — > D* Iv and hadronic modes B° — > 



£>(*)- 



and D* p + for this purpose. The overall effi- 



ciency of our flavor tagging is 99.7% which is consistent 
with the MC expectation. 

Since we know the flavors of both B mesons in this 
case, we can observe the time evolution of neutral B- 
meson pairs with opposite flavor (OF) or same flavor 
(SF), which is given by: 



■PoF(Ai) oc 1 
V SF (At) cx 1 



(1 
(1 



and the OF-SF asymmetry, 
Pof — Psf 



A ■ = 



Vqf + V S i 



2w) cos(AmdAt), 
2w) cos(ATO,jAi), 



1 - 2w) cos(Am d At), 



where w is the wrong tag fraction. We thus obtain the 
value of w directly from the data by measuring the am- 
plitude of the OF-SF asymmetry. 

We obtain the wrong tag fraction by fitting the At 
distribution of the SF and OF events, with Arrid fixed at 
the world average value of 0.472 ps _1 Q. The procedure 
to form the probability density function (pdf ) for the fit is 
quite similar to that adopted for the maximum likelihood 
analysis of CP eigenstates, which is described in the next 
section. 

The resolution function for signal events, which models 
how the true distribution is smeared by the finite vertex 
resolution, is constructed by fitting the proper-time dis- 
tributions without discriminating between the OF and 
SF events and with the lifetime fixed to the world average 
value. In the fit we use the background fraction estimated 
for each region of r, and the proper-time distribution for 
background obtained using events outside the signal re- 
gion. For hadronic modes the sideband regions in 
and AE are used. For semileptonic decays the upper 
sideband in Mdifj is used for the fake D* backgrounds. 
Uncorrelated backgrounds are modeled with the events 
that are found in the signal region after inverting the mo- 
mentum of the lepton. Semileptonic decays D*X£v are 
treated as signal events since they approximately obey 
the same OF-SF asymmetry. 

Figure [ll] shows the measured OF-SF asymmetries as 
a function of At for tagged D* T l^v events for the six 
regions of r. The curves in the figure are obtained by the 
fit. The background is not subtracted in the plots. For 
hadronic modes the fits to OF and SF events are similar 
to those in the semileptonic case. 

We also fit signal MC samples to examine the differ- 
ence between the generated and reconstructed values. We 
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FIG. 11. Measured asymmetries between the OF events 
and the SF events (OF-SF asymmetries) for six regions of r 
obtained for the control sample 5° — > D*iv. The definition 
of the OF and SF events is given in the text. The background 
is not subtracted in the asymmetry plots. Fit curves are also 
shown. 
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TABLE IV. The event fractions (e;) and wrong tag frac- 
tions (wi) for each r interval. The errors include both sta- 
tistical and systematic uncertainties. The event fractions are 
obtained from the J/ipK° s simulation. 



I 



wi(D*£v) ^(hadronic) to; (combined) 



-u.on 

-0.011 

0.149 0.35lig;gr7 



1 0.000 - 0.250 0.405 0.463 

2 0.250 - 0.500 

3 0.500 - 0.625 

4 0.625 - 0.750 

5 0.750 - 0.875 

6 0.875 - 1.000 



0.081 0.254 
0.099 0.169 
0.123 0.107 



+0.021 
-0.020 
+0.019 



+0.015 



0.140 0.0411°;°^ 
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n 24 n+0.021 
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u - 1 '"-0.017 

h0.O22 
-0.014 

041 + 011 
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0.110^ 



apply small corrections to u>z , that correspond to the dif- 
ference. For hadronic modes the corrections range from 
0.003 to 0.03 depending on the r region. For semilep- 
tonic decays the difference is consistent with zero within 
statistical errors, and we apply no correction 

To combine the results from scmilcptonic and hadronic 
decays, we calculate the weighted average and its er- 
ror. We conservatively treat the difference between the 
weighted average and each measurement as an additional 
systematic error, and add this difference in quadrature 
with the error. The systematic errors for the semilep- 
tonic mode are dominated by the uncertainties on the 
background fractions and are compara ble to the statis 
tical errors 



HID 



As explained in Section 
ground estimation relies little on MC information since 
we use control samples whenever possible. One impor- 
tant exception is the C distribution of the D**iv back- 
ground; we use several D** components and add them 
with fixed fractions using MC. Since these fractions are 
poorly known experimentally, we conservatively assume 
that each component dominates the D**tv distribution 
and repeat the fit procedure to obtain the systematic er- 
ror. For the hadronic modes, the main contribution to 
the systematic error comes from the uncertainty of the 
fit bias obtained from the MC simulation, but the sta- 
tistical errors dominate. The event fractions and wrong 
tag fractions are summarized in Table |y| The total ef- 
fective efficiency obtained by summing over the r regions 
is calculated to be 



i 



- 2wi) 2 = (27.0 ± 0.8(stat)t°J(syst))%, 



where ti is the event fraction in each of the six regions. 

Our simulation indicates that events with high- 
momentum leptons dominate the highest r region and 
provide the cleanest tagging information. Events with 
charged kaons have lower r, but are more numerous, and 
thus provide the largest contribution to the effective tag- 
ging efficiency. The effective efficiency using each cate- 
gory alone is examined with MC. We obtain 13% for the 
lepton category, 19% for the kaon and lambda categories 
combined, and 4% for the slow pion category. Note that 



the sum of these values exceeds e c s (29.6% in MC) since 
an event contains tracks in different categories. 

We check for possible biases in the flavor tagging by 
measuring the effective tagging efficiency for the B° and 
B° control samples separately, and for the q = +1 and —1 
samples separately. We find no statistically significant 
difference. 



V. MAXIMUM LIKELIHOOD FIT 

We determine sin20i by performing an unbinned 
maximum-likelihood fit of a CP violating probability 
density function (pdf) to the observed At distributions. 
These pdf's come from the theoretical distributions di- 
luted and smeared by the detector response. For modes 
other than J/tpK*° the pdf expected for the signal is 

e -\At\/r B „ 

V ais {At,q,wi,£f) = — V C p(At,q,wi,£f), (10) 



where 

Vcp{At,q,wi,if) 



2t b o 



1 — £,fq(l — 2wi) sin 2<fii sin(AmrfAt). 

(11) 



In order to take into account the effect of finite ver- 
tex resolution on the At distribution, this pdf is con- 
the back- volved with a resolution function, i? S i g (Ai). Our vertex 



reconstruction method is explained in Section V A . The 
parametrization and extraction of R sl g(At) are described 
in Section V B . We also incorporate the effect of back- 
ground that dilutes the significance of CP violation in 
the time distribution of Eq. ([h]). The At distribution 
for background events, Pbkg(At), is constructed in a sim- 
ilar way to the signal distribution and is described in 



detail in Section V C 



As a result, we adopt the following At distribution 
function for each event: 

P(At i ;sin2^ 1 ) 

= U g J V sig (At',q, wi, ^Rsi^AU - At')dAt' 



+ (1 - / sig )Pbkg(Ai l 



(12) 



where / s i g is the probability that the event is signal, being 
calculated for each candidate from p™ s for J/ipKl and 
a combination of AE and M^ c for other modes. The 
only free parameter is sin 2<f>i , which is determined by 
maximizing the likelihood function 



L = Y[P(AU;shx2(pi), 



(13) 



where the product is over all candidates. We perform 
a blind analysis: The fitting algorithms were developed 
and finalized without using the flavor information q. 

In the following we explain the details of i? s j g (Ai), 
Pbkg(Ai) and / sig in turn. The likelihood for J/ipK*°{^ 
KgTr ) candidates is described separately in Section VE. 
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A. Vertex Reconstruction 

The decay vertices for the CP side that include a J/ip 
candidate are reconstructed using leptons from the J/ip 
and a constraint on the B decay point. The B decay 
point is constrained by the measured profile of the in- 
teraction point (IP profile) convolved with the finite B 
flight length in the plane perpendicular to the z axis (the 
r-(f> plane). The IP profile is represented by a three- 
dimensional Gaussian distribution. The standard devi- 
ation of each Gaussian is determined using pre-selected 
BB candidates on a run-by-run basis, while the mean is 
evaluated in finer subdivisions. The typical size of the IP 
profile is 100 ^m in x, 5 /zm in y and 3 mm in z. Since the 
size in the y direction is too small to be measured from 
the vertex distribution, it is taken from special measure- 
ments by the KEKB accelerator group. For leptons, we 
require that there are sufficient SVD hits associated with 
a CDC track by a Kalman filter technique; i.e. both z 
and r-cf) hits in at least one layer and at least one ad- 
ditional layer with a z hit. In order to remove events 
with mis-reconstructed tracks, wc require that the re- 
duced x 2 (x 2 / n , n = number of degrees of freedom) of 
the vertex be less than 20. The vertex reconstruction 
efficiency is measured to be 95% with B ± — ► J/ipK^ 
and B° — > J/ipK*°(-^> K ± ir T ) events. This is consistent 
with the expectation from the SVD acceptance and clus- 
ter matching efficiency. The resolution estimated by MC 
is typically 75/zm (rms). 

For B — > r] c Kg candidates, the method is basically the 
same as for J/ip, replacing ( + (r with K + K~ in case of 
Vc -» K+K-ir and with K ± n^ for Vc -» K^K***. 
Although the resolution in these cases is worse than for 
candidates with a J/ip vertex, it is still better than the 
tag-side vertex resolution. 

The algorithm for tag-side vertex reconstruction is cho- 
sen to minimize the effect of long-lived particles, sec- 
ondary vertices from charmed hadrons and a small frac- 
tion of poorly reconstructed tracks. From all the charged 
tracks except those used to reconstruct the CP side, we 
select tracks that have associated SVD hits in the same 
way as for the CP side. We also require that the im- 
pact parameter with respect to the CP-side vertex be 
less than 0.5 mm in the r-<j> plane, less than 1.8 mm 
in z, and the vertex error in z be less than 0.5 mm. 
Tracks are removed if they form a Kg candidate satis- 
fying | M K o - M n+n -\ < 15 MeV/c 2 . Tracks satisfying 
these criteria are used to reconstruct the tag-side vertex 
where the IP constraint is also applied. If the reduced x 2 
of the vertex is good, we accept this vertex. Otherwise 
we remove the track that gives the largest contribution to 
the \ 2 and repeat the vertex reconstruction. If the track 
to be removed is a lepton with p^ ms > 1.1 GeV/c, how- 
ever, we keep the lepton and remove the track with the 
second worst x 2 ■ This trimming procedure is repeated 



until we obtain a good reduced x 2 ■ The reconstruction 
efficiency was measured to be 93% for B ± — ► J/ipK ± and 
B° — > J/ipK*°(^ K + ir~) candidates, consistent with 
the MC expectation (91%). The resolution estimated 
from the simulation is typically 140/xm (rms). 

B. Signal Resolution Function 

The resolution function i? s ig(Ai — At') is parametrized 
by the sum of two Gaussians: 

R sig (At - At') = (1 - / tai i)G(Ai - At'; hau a At ) 

+ / tail G(At-At'; M t A a j 1 , C 7if), (14) 

where G(x; fx, a) is a Gaussian distribution in x with 
mean fj, and rms a. The parameter / ta ;i describes the 
fraction of the tail of the resolution function, and a At, 
^At'i Ma* and y^^t are the proper-time difference resolu- 
tions and the mean value shifts of the proper-time dif- 
ference for the main part and the tail of the resolution 
function, respectively. The value of / ta ii is determined 
to be 0.03 ± 0.02 from the lifetime analysis of hadronic 
samples using the same resolution function. 

The proper-time difference resolutions aAt and a^f 
are calculated on an event-by-event basis taking into 
account the error uk in the kinematic approximation 

/ /T tai ' 

^t=\j(^r+^r) 2 - 

We measure a K = 0.287 ± 0.004 ps and a^ n = 0.32 ± 
0.19 ps using the MC simulation. These parameters are 
independent of the detector performance. 

The parameters oaz and are calculated from the 
event-by-event vertex errors of the two B mesons, a^' p 
and er' ag , which are computed from the track helix errors 
in the vertex fit. We use 

2 _ q2 („CP\2 , ( q2 , q2 \(„ta%\2 
a Az — D dctl°z ) "+" Wdot ~r J charmJl CT z ) ! 

(a^f = {S$S ) V? P ) 2 + {(5^) 2 + (S^iJ 2 }^) 2 , 

where S'charm and S^arm are scaling factors to account 
for the degradation of the vertex resolution on the tag- 
side due to contamination from charm daughters, and 
Sdot and S^ft are global scaling factors that account 
for systematic uncertainties in the vertex errors a^ p 
and CT* ag . We determine S'charm = 0.58 ± 0.01 and 
S^him = 2.16 ± 0.10 using the MC simulation. The val- 
ues of Sdet and are measured from the data as they 
depend on the detector performance. We determine Sdet 
using a D° — > K + ir~ control sample. The production 
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point of the D° is obtained from the primary tracks in 
the same hemisphere as the D° candidate using the IP 
constraint. The distance between the D° decay vertex 
and the production vertex in the z direction is fit with 
the same resolution function and the known D° lifetime 
to obtain Sdet- We measure S^et = 0.88 ± 0.01 from 
the data and S^ct = 1-05 ± 0.01 from a MC simulation 
of the D° sample. Finding S dct = 1-035 ± 0.003 for a 
B J/ipK MC sample, we use Sdet = (0.88 ± 0.01) x 
(1.035 ± 0.003) /(1.05 ± 0.01) = 0.86 ± 0.01 for the data. 
We determine to be 3.51 ± 0.88 from the lifetime 
analysis of the flavor-specific hadronic samples. 

For B — > r) c Kg decays, we introduce an additional 
scale factor to account for the difference between the 
r/ c and J/ip decays. Using the MC, we determine the 
additional scale factor for r\ c to be 1.07 ± 0.02 for both 
r) c -► K+K-ir° and rj c -> K^K^. 

A small fraction of events have a large reduced \ 2 ■ We 
have found that the vertex error computed from track 
helix errors in the vertex fit underestimates the vertex 
resolution and the vertex with larger \ 2 has worse res- 
olution. In order to take into account this effect, we 
introduce effective vertex resolutions and <r* ag when 
X 2 /n is greater than 3: 

fc CP ) 2 = [l + Q CP{(x»cP-3}]k CP ) 2 , 

(af g ) 2 = [1 + ata g {(x 2 /«)t ag - 3}](a^) 2 , (15) 

where (x 2 / n )cp and (xV n )tag are the reduced \ 2 01 
the vertex fits for the CP and tagging B decay vertices, 
respectively. The coefficients acp — 1.02 ± 0.03 and 
"tag = 1-64 ± 0.05 are determined by a MC study of 
B° -> Jj^K%. 

As mentioned above, the offsets fiAt and /i^J 1 originate 
from the mean shifts of the Az measurements /ja z and 
Ma*) respectively: 

HAz 

c(3j' 



fJ-At 



tail 

Ma* 



. .tail 
MAz 

c(3j ' 



The mean value shifts, fiAz and fJ,^, are caused by con- 
tamination from charm daughters in the vertex recon- 
struction on the tag-side and are correlated with cr* ag : 



tag 

z > 



MAz(c z as ) = jito + a^a 

..tail/' tag\ ,,tail , „tail tag 
MAz \°z ) = M0 + a n a z 



The values for no and /Xg al1 are determined from hadronic 
samples to be ^ = (—21.4 ± 3.7) ^m and /ip al1 = 
(151 ± 128) /im, while and ajf 1 ' are derived from 
MC simulation where we obtain = —0.10 ± 0.01 and 
a* 8 " = -1.42 ±0.17. 

Figure [12 (a) shows the Ai rec — At gen distribution for 
the MC J/ipKg(ir + Tr~) candidates along with the res- 
olution function, where At rcc and At gcn are the recon- 
structed and true proper-time differences, respectively. 



The resolution function is obtained by summing event- 
by-event resolution functions. The distribution is well- 
represented by the resolution function. The average res- 
olution function obtained from the J/ipKg(ir + Tr~) data 
is shown in Fig. [l^ (b), which is represented by the sum 
of two Gaussian distributions with the following param- 
eters: A*main= -0.24 pS, Mtail= 0.18 pS, CT main = 1.49 pS, 
Ctail= 3.85 ps. 



C. Background Shape 

The background likelihood function is defined in a sim- 
ilar way to the signal function, 



fbkg(At) 



Pbkg(At') • i? bk g(Ai - At')dAt'. (16) 



Although i?bkg(Ai — At') is treated as a resolution func- 
tion for the background, it does not need to be an exact 
model of the vertex resolution. It is more important that 
■Pbkg represents the proper-time distribution of the whole 
background sample with sufficient precision. 
The pdf for background events is expressed as 

— |A*l/ r bkg 

V hig (At)=f T — + (1 - f T )S(At'), (17) 

^Tbkg 

where f T is the fraction of the background component 
with an effective lifetime of ibkg, and 5 is the Dirac 
delta function. We assume no asymmetry in the back- 
ground At distribution. We find f T to be small using 
background-dominated regions in the AE versus M bc 
plane of J/tpKg and J/ipK + candidates. We thus use 
■pbkg(Ai) = 5 (At) for all fcp modes except for J/ipK®. 

The background in the J/tpK® mode is dominated by 
B — > J/ipX decays, including CP eigenstates that have 
to be treated differently from non-CP states. The Pbkg 
for the J /tpKl mode is determined by a MC simula- 
tion study separately for each background component: 
J/iPK*°(KItt°), = -1 CP modes (J/ipK° s ), = +1 
CP modes (^(2S)K° L , XdK Q L and J/ijm°), the other B°, 
B ± decays and combinatorial background. For the CP- 
mode backgrounds, we use the signal pdf given in Eq. (|ll| ) 
with the appropriate £/ values. For the J/t/jK*°(Kl7T°) 
mode, which is a mixture of £/ = — 1 (about 81%) and 
£/ = +1 (about 19%) states j33[, we use a net CP eigen- 
value of ^ K ,o = -0.62 ± 0.11. 

Accordingly, we obtain the background pdf for J/ipK^ 



Pbkg (At) 



x i fb g B° 



2t b o 

+fb g ^K'°'Pcp(At, q, wi,£^ K *o) 
+fb g cp odd 'Pcp(At,q, wi, -1) 

+fb K CP cvcn 'Pcp(At, q, wi, +1) >+/b g B± 



(18) 

e -|At|/r bgB ± 

2r b gB± 
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e ~|At|/r bg Cmb 

"/bgCmb i /rCmb ^ ^ Q ~ /rCmb)^(At) ^ , 



2r b , 



Cmb 




At rec- At gen (PS) 



10 



>40 

1 

< 1 



10 



■ (b) 


/-X — main + tail 


, , , , i , , 


/ \ ---- tail 

, , i i i ■ ■ i i i i i 



-10 -5 5 

At'-At (ps) 
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FIG. 12. (a) The At [cc - At gen distributions for the MC 
candidates with the resolution function and 
(b) the average resolution function from J/iPK s (tv + tv~) data. 
The vertical scales are arbitrary. 



where fbgfK- , fb s CP odd , /b g CP ovonJ fb g B°, fbgB±, and 
/bgCmb are the fractions of background components from 
J/iPK* Q {K° L ir°), £f = -1 CP-modes, = +1 CP- 
modes, the remaining B°, B ± , and combinatorial re- 
spectively (fb g ^K*o + /bgCP od d + /b g CP ovo „ + fb S B° + 
/bgB± + /bgCmb = !)■ The fraction of each back- 
ground component is a function of p™ ls . The fraction, 



/bgCmb, is calculated as described in Section [III Q . The 
combinatorial background includes a prompt component 
with the fraction of (1 — / T cmb) where / T Cmb = 0.26 
± 0.08. The lifetime distribution of the combinatorial 
background is obtained from e-/i combinations with in- 
variant masses in the J/ip region that satisfy our selec- 
tion criteria. The effective lifetime is determined to be 
T bgCmb = 1-03 ± 0.12 ps, also from the e-fi control sam- 
ple. A MC study shows that the effective lifetime for 
background from B ± , Tb g s±, is shorter than the B ± life- 
time due to the contamination of charged tracks from 
fcp (mostly 7r ± from J/ipK* ± (Klir ± )) into the tag-side 
vertex. The value of r bgB ± is determined from the MC 
simulation to be (1.49 ± 0.04) ps. The same MC study 
shows that the effective lifetime for B° backgrounds is 
consistent with the nominal B° lifetime. Thus we use 
the nominal B° lifetime in our fit. 

For the J/tpX background in the J/tpK^ mode, we 
use the signal resolution function R slg to model the back- 
ground since both the CP- and tag-side vertices are re- 
constructed with similar combinations of tracks for these 
backgrounds. 

For the combinatorial background, we use 



P bkg (A<'- At)= (1-/, 



tail.bk; 



;)G(Ai' 



<\ f . bkg bkgN 

' Ma* > "At ) 



+/taii,bkgG(At'-Af;^t bkg , 



tail, bkg 
r At ' 



). (19) 



where /taikbkg, P^t 



j tail.bke 

and fj, At 



are constants deter- 



mined from data. The resolutions and a 



calculated on an event-by-event basis as 



tail. bkg 
7 At 



are 



bkg 



tail, bkg 
^ At 

/ tail, bkg \ 2 
\ a Az ) 



a 



bkg 



c/?7 

tail. bkg 



a 



Az 



c/37 

(s'btn^y + ia^) 2 } 



where a^ p and cr' ag are calculated as shown in Eq. (|lT_ 

We use different values of ^ k t s and /ti^t ' 8 for the fi- 
nite lifetime component and the zero-lifetime component, 
since they come from different types of events. The back- 
ground shape parameters for all modes except J/ipK^ 
are obtained from events in the background-dominated 
regions of AE versus M^ c . For J/ipK®, we use events 
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with e-fi pairs to determine the properties of fake J/ij) 
candidates, as discussed in Section [II C . The parame- 
ters used in the fit are summarized in Table |v|. 

TABLE V. Background shape parameters for the combi- 
natorial background. 



parameters 


ccK°(K*°) 




obkg 
°dct 


1.08 ±0.06 


0.88 ± 0.07 


otail 
°bkg 


3.31 ±0.28 


3.94 ± 1.14 


J^tailjbkg 


0.14 ±0.04 


0.05 ± 0.02 


lifetime component 






MAt g (P S ) 


N/A 


-0.33 ±0.10 




N/A 


1.86 ± 1.28 


prompt component 






Ma?(P s ) 


-0.05 ± 0.04 


-0.22 ± 0.15 




-0.12 ±0.26 


-5.00 ± 1.10 



D. Signal probability 

The signal probability, / s i g , is calculated as a function 
of AE and Mb c for each event. It is given by 

F SIG (AE,M hc ) 



f sig (AE,M hc 



F BG (AE, M bc ) + F SIG (AE, M bc 



(20) 



and 



where Fsig(A_B, M\, c ) is the signal function 
Fbq(AE, Mbc) is the background function. 

In the case of B° — > J/ipKs each distribution of AE 
and Mb c is well modeled by a Gaussian function. For 
the background, we use a linear function for AE and the 
ARGUS parametrization |28| for Mb c : 

Fsig(AE, Mbc) 

= a ■ G(AE\ (iaeiVae) ■ G(M hc ; ^ Mbc , er Mbc ), 
Fbg (AE, Mbc) 

= b ■ (1 + C • AE) ■ M bc y/1 - (Mbc/^bcam) 2 

x exp(n • [1 - (M bc /£ b cam) 2 ]) (21) 

where a and b are normalization factors consistent with 
the overall signal-to-background ratios obtained from the 
ht to the Afbc distribution in the AE signal region. The 
values u Af bc ; ft > a M^ E ^ MM bc > c and n are determined 
from a fit to the data. 

The AE and M bc distributions for B° -> 
V>(2S)tf s (V>(2S) -» ^ + n, S° ^(2S)# S (^(2S) 
J/^tt+tt-), and B° -► rycXg (r? G -> i^i^TrT), are de- 
termined using the same procedure as that for B° — > 
J/ipKs- For B° ->■ Xci-^s, the AS and M bc distribu- 
tions are determined from MC simulation because the 
data sample for this mode is too small to estimate the 
parameters reliably. 



The treatment for modes that include ir° mesons such 
as B° -> J/iPK s (K° s -> 7r°7r°) and B° -> r) c K° s (r/ c -> 
K + K~7T°), is different. While the fit function for the 
M bc distribution remains the same, the AE distributions 
are better represented by the Crystal Ball function [M : 



A eX Pl 

x exp(-q 2 /2) 
A 



aaAE 



for AE > fjL^E 
for AS < hae — okjae- 



All the parameters for these fits were determined from 
MC simulation because the number of events for these 
modes in data is too small. The integrated background 
fractions in the signal region are listed in Table VI. 

For the B° — > J/ipK { l fit, we define the signal proba- 
bility as a function of p™ s , as described in Section III C . 



TABLE VI. Summary of the numbers of candidates and 
background fraction in the signal region for each mode. The 
values are obtained for events that have successful vertex re- 
construction and flavor tagging. 



Decay mode 



B u 
B° 
B° 
B° 
B° 
B° 
B° 
B° 
B° 



K 



JJWs, 

J/i>K° s , K° s -+ ttV 

ip(2S)K° s , 4>{2S) -» J/i/>tt h 

Vc K° s , nc -» K° S K±^ 
n c K%, rj c -> K + K~n° 

■ J/i>K° L 

J/i>K*°, K*° -» K° s n° 



Events bkg. fraction 



387 
57 
33 
32 
17 
35 
17 

523 
36 



0.038 ± 0.010 
0.272 ± 0.054 
0.038 ± 0.028 
0.078 ± 0.027 
0.144 ± 0.056 
0.242 ± 0.045 
0.560 ± 0.164 
0.379 ± 0.048 
0.163 ± 0.054 



E. Likelihood for J/ipK*°(-^> K%tz°) 

For the B° — ► J/tpK*° fit, the signal pdf we use is: 

e -\At\/r B o 

r sig (At,e tI ,q,wi,£,f) = — x 

Zt b o 

[(1 - /odd)g(l + cos 2 9 tI )V C p(At, q,w l7 +l) 

+/odd|(l - cos 2 9 tl )Vcp(At,q, Wl ,-l)}, (22) 

where / 0( jd is the fraction of = —1 decays in the 
B° -> J/iJjK* (K* -> K°ir°) mode determined from 
a full angular analysis to be 0. 19 ± 0.04 ± 0.04 ||. Here 
6* tr is defined in the transversity basis |35| as the angle 
between the positive J/tp decay lepton direction and the 
axis normal to the K*° decay plane in the J/ip rest frame. 
Vcp{At,q,u>i,£f) is defined in Eq. (JTT|) . The signal res- 
olution function is identical to that used for the other 
modes. For the background shape, we also use Eq. ( fl9| ) 
for i? b kg except for the J / ipX background where we use 
i?si g in the same way as for the J/ipK® fit. 
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We use the following background pdf: 



TABLE VII. Summary of sin20i fit results. Only statis- 
tical errors are shown. 



Pbkg(At) = /bgFA 



,-\At\/r B 



,-\At\/TB Nn 



/bgCmb<5(At), 



/bgNR- 



(23) 



where /bgFA, jbgNR, and /b g Cmb are the fractions of back- 
ground components from feed-across from other J/ipK* 
modes, non-resonant B° — > J/ipK s TT° decays and com- 
binatorial background. The fractions of feed-across and 
non-resonant decays are determined from the MC simu- 
lation and from K* mass sideband events, respectively, 
and are functions of M\ x . The fraction of combinatorial 
background is determined in the same way as for the K s 
modes. The effective lifetimes of the feed-across and non- 
resonant decay backgrounds, tb fa and tb nr , are fixed to 
the B° lifetime in the fit. 

Finally, the determination of / s i g follows the method 
for other modes that include ir° mesons. The AE dis- 
tribution is modeled by a Crystal Ball function. We 
consider contributions from the feed-across from other 
B a — > J/ipK*° modes as well as from the non-resonant 
B° — > J/t/jKgTT mode, which make a peak in the sig- 
nal region. These background fractions are determined 
from the MC simulation and K*° mass sideband data, 
respectively. 



VI. FIT RESULTS 

The likelihood fit is applied to the 1137 candidates 
where the vertex reconstruction and flavor tagging have 
been successful. We obtain 

sin 2<?(>i = 0.99±0.14(stat). 

The observed CP violation is large. Figure [l3| shows 
the At distributions together with the results from the 
fit. Indeed, the broken CP symmetry is visually appar- 
ent from the difference between the number of events for 
q£f = +1 and q£f = —1 at each At bin, despite the dilu- 
tion from the vertex resolution, background events and 
incorrect flavor tagging. 

We examined the value of sin 20i in various sub- 
samples. Applying the likelihood fit to (cc)K s (£/ = — 1) 
and J/ipK® (£/ = +1) separately, we obtain sin20i = 
0.84±0.17(stat) and 1.31 ± 0.23(stat), respectively. Fig- 
ure [l4| shows the log-likelihood values as a function of 
sin 20! for CP-odd, CP-even, and all decay modes. A 
more detailed breakdown along wi th se parate results for 
q = +1 and —1 is given in Table VII. We find no sys- 
tematic trends beyond statistical fluctuations. 

Figure [l^ shows the asymmetry, sin 20i • sin(Am,jAt), 
obtained in each At bin for (a) all modes, (b) CP-odd 
modes and (c) CP-even modes. The unbinned maximum 
likelihood fit is performed separately for events in each 



Sample 


Events 


sin 20i 


/tag = B u (q = +1) 


560 


0.84 ±0.21 


/tag — o \q— rj 


^77 


1.11 ±0.17 


(cc)K» s 


578 


0.84 ±0.17 




387 


0.81 ±0.20 


(cc)K° s except J/%})K s (it + k-) 


191 


1.00 ± 0.40 


J/^K L 


523 


1.31 ±0.23 


J/^K*°(K° s -k ) |36| 


36 


0.97 ± 1.40 


All 


1137 


0.99 ±0.14 



At bin. The value of sin 20i and its error are multiplied 
by the average value of sin( Arrid At) in each At bin of the 
plot. The points are plotted at the average At of each 
bin. 



0.20 



< 

T3 



0.10- 



0.00 




FIG. 13. At distributions for the events with q£f — ±1 
(solid points) and q£f = —1 (open points). The results of the 
global fit (with sin 20i = 0.99) are shown as solid and dashed 
curves, respectively. 



We also checked the values of sin20i in the different 
r rang es of the flavor tagging. The results are listed in 
No systematic variation is seen. 



Table VIII . No systematic variation is seen. Finally, 
we subdivided the CP sample into three data taking 
periods: in 2000, from January to April 2001 and the 
rest. The sin20i values we obtain are 0.84 ± 0.29(stat), 
1.33±0.30(stat) and 0.94±0.20(stat), respectively. Again 
the results are consistent within the statistical fluctua- 
tions. 

TABLE VIII. The r dependence of sin 20i fit result. Only 
statistical errors are shown. 



r region 


0.0-0.5 


0.5-0.75 


0.75-0.875 


0.875-1.0 


Events 
sin 20i 


613 

o.6o ±g : » 


239 
0.40 t° H 


119 

I r fi +0.25 
l.OO _ .29 


166 
1 05 + - 15 

l - uo -0.18 
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CM 

I 




sin2<)) 1 

FIG. 14. Values of — 21n(L/L max ) vs. sin20i for the 
£y = —1 (dashed) and +1 (dotted) modes separately and for 
both modes combined (solid). 



A. Systematic Errors 

The sources of systematic error we consider are listed 
in Table IX pTJ. Adding all the systematic errors in 



quadrature, we obtain 

sin 2<?(>i = 0.99 ± 0.14(stat) ± 0.06(syst). 
Below we explain each item in order. 

TABLE IX. List of systematic errors on sin 20i . 



source 


+error 


—error 


vertex reconstruction 


+0.040 


-0.040 


resolution function 


+0.022 


-0.032 


wrong tag fraction 


+0.022 


-0.025 


physics (t s o, Ami, J/tpK*°) 


+0.007 


-0.004 


background fraction (except for J/tpK^) +0.003 


-0.004 


background fraction (J/tpKl) 


+0.020 


-0.020 


background shape 


+0.001 


-0.001 


total 


+0.06 


-0.06 





1 - 





-1 - 




-8 



(a) Combined 



(b) (cc)K s ^ f =-V t 



(d) Non-CP sample 



• • •- 



4 
At (ps) 



1. Vertex reconstruction 

The largest contribution comes from vertex reconstruc- 
tion. We searched for possible biases by using two differ- 
ent vertexing algorithms and changing the track selection 
criteria for the tag-side vertex. In the alternative vertex- 
ing algorithm, we first obtain a seed vertex using tracks 
of good quality: an impact parameter from IP in r-<f> di- 
rection is smaller than 2.5 times the r-<p vertex error; the 
vertex error in z is less than 0.5 mm; and the cms mo- 
menta are larger than 0.3 GeV/c. We then repeat the 



FIG. 15. (a) The asymmetry obtained from separate fits 
to each At bin for the full data sample; the curve is the re- 
sult of the global fit. The corresponding plots for the (b) 
(cc)K° s = -1), (c) J/i>Kl {i f = +1) and (d) non-CP B° 
control samples are also shown. The curves are the results 
of the combined fit applied separately to the individual data 
samples. 
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vertex fit using tracks within 3 (4) a in z from the seed 
vertex for the cms track momentum less (larger) than 1 
GeV/c, where a is the error of the seed vertex in z. We 
also estimated the effects of the vertex resolution tails 
using samples with small At (\At\ < 6 ps) and tighter- 
vertex quality cuts. 



2. Resolution function 

We estimate the contribution due to the uncertainty in 
the resolu tion function by varying its parameters (given 
in Section |VB|) by ±la. 



6. Background fraction for J /ijjK^ 



3. Wrong tag fraction 

Systematic errors due to uncertainties in the wrong 
tag fractions given in Section IV B (Tabic IV) are stud- 
ied by varying wrong tag fraction individually for each r 
region. We added the contributions from each variation 
in quadrature. 



4- Physics parameters 

The B meson lifetime and mixing parameter are fixed 
to the world average values in our fit; i.e. tb = 
(1.548 ±0.032) ps and Am d = (0.472 ± 0.017) ps" 1 . We 
estimate the systematic error by repeating the fit varying 
these parameters by their errors. Another physics-related 
uncertainty is the CP eigenvalue of J/ipK*° (£j/^k'°) 
measured from the angular distribution of the decay 
daughters j33). This systematic uncertainty is deter- 
mined from the ±lcr uncertainty in the measurement. 



As described in Section III C[ , the background fraction 
for the J/ipK^ sample is obtained from a fit to p™ s dis- 
tribution and is given in Table 0. In this fit, the sum 
of components is automatically constrained to the total 
number of events in the signal region. Thus, the signal 
yield and the size of other backgrounds are strongly anti- 
correlated. To determine the systematic error on sin 2</>i 
that comes from the uncertainty of the background, we 
need to take this anti-correlation into account. To this 
end, we repeat the fit to the p™ s distribution with the 
background fractions as free parameters but with the sig- 
nal yield fixed +ler or —la away from the central value 
obtained in the nominal fit. The resultant background 
yields are used to repeat the procedure to obtain sin 20i . 
We regard the difference between the thus obtained value 
and our nominal sin 20i value as the systematic error. We 
also check the systematic error due to the uncertainty in 
the CP content of the background. We repeated the fit 
varying parameters to determine the various background 
fractions. Since these parameters are obtained from the 
MC simulation, we estimate the systematic error by con- 
servatively changing each parameter by ±2tr and adding 
the resulting changes in quadrature. 



7. Background shape 

The parameters that determine 'Pbkg(At) and i?bkg, 
given in Section V C , are varied within their errors and 
fits are repeated. 



B. Cross checks 



5. Background fraction except for J/ipK l l 



The background fraction in our pdf, 1 — / S j g , is calcu- 
lated from the signal and background distribu tion func- 
tions of AE and Mb c as described in Section VP . The 
distribution functions of AE and M^ c 



are determined 
from data or the MC simulation depending on the decay 
mode. To estimate the systematic errors associated with 
the choice of parameterization, we varied the parameters 
obtained from the MC simulation by ±2cr and the pa- 
rameters obtained from the data by ±1ct. The likelihood 
fit was repeated. A wider range of uncertainty was con- 
servatively chosen for parameters obtained from the MC 
simulation to take into account the possible difference 
between the MC simulation and data. We also estimated 
the systematic errors for the integrated background frac- 
tions, listed in Table by varying these parameters by 
±1(7. We added the results of these calculations for each 
decay mode in quadrature. 



We performed several cross checks: The fitting proce- 
dure was examined using MC samples based on our likeli- 
hood functions (toy MC samples) . We also measured the 
B meson lifetime using the same vertex reconstruction 
method. In addition we tested non-CP control samples. 
These cross checks are described below. 



1. Ensemble test 

A thousand toy MC samples, each containing 1137 
events, are generated based on our likelihood function 
to check the fitting procedure. Figure |l6| shows the dis- 
tribution of residuals (sin20i(fit) — sin 20i (input)), pulls 
(residuals divided by fit errors), and the positive and neg- 
ative errors on sin20i returned from the fits. All the toy 
MC samples have an input value of sin20i = 0.99. The 
center of the residual in Fig. [l6| (a) is consistent with 
zero, and the standard deviation of the pull distribution 
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in Fig. |l6| (b) is consistent with unity. Therefore the 
global fit returns the input sin 2(f>i value and a reasonable 
error. Figure [l6| (c) and (d) also show that the positive 
and negative errors obtained from the fit are consistent 
with expectations. 




positive error negative error 



FIG. 16. The result of our toy Monte Carlo test 
of our fit: (a) the distribution of residuals ( sin 2(/>i (fit) 
— sin 20i (input) ), (b) pull distribution ( (residual) /(error of 
fit) ), and (c) the positive and (d) negative errors of sin20i 
returned from the fits. The vertical lines in figures (c) and 
(d) indicate the errors obtained from the fit to data. 

We also generate toy MC samples for J/ipKg and 
J /ipK^. The average negative errors are 0.19 and 0.28 for 
J /ip Kg and J/ipK^, respectively. Our measurement for 
J/ipKg (0.20) is in good agreement with the expectation, 
while the result for for J/iftK^ (0.23) is smaller. We ob- 
tain the probability of obtaining smaller errors than this 
measurement to be 1.4% for J/ipK^, which is within a 
possible fluctuation. 

2. B° lifetime 

The B° lifetime has been measured with the same data 
sample. We apply the same vertex reconstruction algo- 
rithm for fully reconstructed B° decays as for the CP 
decays and the tracks on the tag-side. Unbinned max- 
imum likelihood fits are made with an exponential pdf 
convolved with the same At resolution function and back- 
ground pdf as in the fit for CP eigenstates. For the 
combined B° -> J/if)K*°(K + TT-), D^-ir+, D*~ p+ and 
D*~l + v decay modes, the B° lifetime is measured to be 
t b o = 1.547 ± 0.021(stat.) ps. The result is consistent 
with the world average value [[| . 



3. Tests on control samples 

We use control samples of non-CP eigenstates, B° — > 
J/ipK*(K + Tr-), B° -► DW-7T+ B° -> D*-p+ and 
B° — > D*~ £ + v, to check for biases in the analysis. We 
perform the same fit to these control samples as for the 
CP-eigenstate modes. The results, summarized in Ta- 
ble |x|, show no systematic tendency. A combined fit to 
all the modes yields 0.05 ± 0.04, consistent with zero at 
the 1.2er level, as shown in Fig. [l5| (d). 



TABLE X. CP asymmetries for control samples. 





J/4<K* (K + tt-) 


D ( * ) -tt+ D*~£+u 






and D*~ p + 


Events 


816 


5560 10232 


CP asymmetry 


0.01 ±0.14 


0.12 ±0.06 0.01 ±0.05 



We check for a possible bias due to CP asymmetry in 
the background. We fit J/tpK^ candidates in the back- 
ground region (1.0 < p™ ls < 2.0 GeV/c) treating all the 
events as J/tpK^ candidates. Note that the fraction of 
events with definite CP in this region of p™ s is expected 
to be negligible. The result is sin20i = 0.49 ±0.35, 
consistent with zero at the 1.4cr level. 

C. Discussion 

We have performed several statistical analyses of the 
results described in the previous sections. Using a Gaus- 
sian likelihood function based on the statistical and the 
systematic errors, we calculated the confidence intervals 
bounded by the physical region for sin 20i using two 
methods: the Feldman- Cousins Q frequentist approach 
and the Baycsian method with a flat prior pdf. We find 
a lower bound on sin20i of 0.70 at the 95% C.L. in 
both cases. We also estimated the Bayesian lower limit 
using the exact likelihood function, shown in Fig. |l4| , 
and obtained 0.69. We conclude that the likelihood 
function is Gaussian to a good approximation. Com- 
binations of indirect measurements typically constrain 
0.50 < sin 20! < 0.86 in the framework of the SM |9). 
Although our measured value is large, it is consistent 
with the higher range of the SM prediction. We are con- 
tinuing the measurement with much higher statistics in 
order to test the KM ansatz more precisely. 

Finally we comment on the possibility of direct CP 
violation. The signal pdf for a neutral B meson decaying 
into a CP eigenstate (Eqs. |l0| and |ll|) can be expressed 
in a more general form as 

e -|At|/T B o 

f 1±JA|_ _ ^ _ 2wi)MM sin ( Amd At) 
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|A| 2 



■ cos(AmdAt)] 



(24) 



where A is a complex parameter that depends on both 
B°B° mixing and on the amplitudes for B° and B° decay 
to a CP eigenstate ||. The coefficient of sin(AmdAi) is 
given by _4|A| = — £/iroA and is equal to sin2</>i in the 
SM. The presence of the cosine term (|A| =^ 1) would 
indicate direct CP violation. Throughout our study we 
have assumed |A| = 1, as expected in the SM Q. In order 
to test this assumption, we also performed a fit using 
the above expression with A and |A| as free parameters, 
keeping everything else the same. We obtain |A| = 1.09± 
0.14 and A = 0.80 ± 0.19 for the J/ipK%(n + n-) sample 
only, and |A| = 1.03 ± 0.09 and A = 0.99 ± 0.14 for 
all CP modes combined, where the errors are statistical 
only. This result confirms the assumption used in our 
analysis. 



VII. CONCLUSION 

We have measured the CP violation parameter sin 20i 
at the KEKB asymmetric e + e~ collider using a data sam- 
ple of 29.1 fb _1 recorded on the T(4S) resonance with 
the Belle detector. To extract sin20i, we apply a maxi- 
mum likelihood fit to the 1137 candidate B meson decays 
to CP eigenstates. We obtain 

sin 20i = 0.99 ± 0.14(stat) ± 0.06(syst). 

We conclude that there is large CP violation in the neu- 
tral B meson system. A zero value for sin2</>i is ruled 
out by more than six standard deviations. 

We have calculated the confidence intervals bounded 
by the sin20i physical region for the Bayesian and 
Feldman-Cousins frequentist approaches. The lower 
bound is found to be 0.70 at the 95% C.L. in both cases. 
Our result is consistent with the higher range of values 
allowed from the constraints of the Standard Model as 
well as with our previous measurements B . 
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